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Abstract 



The surface phonon dispersion curves of commensurate Xe monolayers on 
Cu(lll) and incommensurate Xe monolayers on Cu(OOl) surfaces have been 
measured using He atom scattering (HAS) time of flight (TOF) spectroscopy. 
The TOF spectra are interpreted by combining quantum scattering calcula- 
tions with the dynamical matrix description of the surface vibrations. Both 
a vertically polarized Einstein-like mode and another, acoustic-like mode of 
dominantly longitudinal character, are identified. The latter mode is charac- 
terized by the presence and absence of the zone center frequency gap in the 
commensurate and incommensurate adlayers, respectively. The microscopic 
description of the TOF spectral intensities is based on the extensive theoret- 
ical studies of the interplay of the phonon dynamics, projectile-surface po- 
tentials, multi-quantum interference and projectile recoil, and their effect on 
the HAS spectra. Both single and multi-quantum spectral features observed 
over a wide range of He atom incident energies and substrate temperatures 
are successfully explained by the theory. 
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I. INTRODUCTION 



In the past, investigations of physisorption potentials for particles on solid substrates as 
well as the potentials describing interadsorbate interactions were mainly motivated by the 
experimental studies of adsorption and phase transitions at surfaces. Several new develop- 
ments have recently stimulated renewed interest. First, the possibility to manipulate and 
displace single atoms and molecules with low-temperature scanning tunneling microscopes 
||T|J^ has not only added to our understanding of adsorbate-substrate interactions but also 
underlined the need for more detailed microscopic data. Second, the recent progress in 
understanding sliding friction |^ calls for more precise information on the potential energy 
surfaces and energy dissipation to the substrate [|]. Although physisorption energies can 
be determined in a rather straightforward fashion using thermal desorption spectroscopy 
(TDS), information on the interparticle potentials can be only extracted indirectly, either 
through precise analyses of the shape of TDS-peaks or by detailed investigations of phase 
transitions within the adsorbed adlayers 0. 

As the interplay between adsorption and interadsorbate potentials determines the low 
energy dynamics of adlayers, studies of physisorbed rare gas monolayers are of special interest 
in the above discussed context because of their allegedly simple vibrational properties. In 
the present paper we describe a determination of the parameters of these potentials on the 
prototype commensurate and incommensurate physisorbed monolayers of Xe on Cu(lll) and 
Cu(OOl), respectively, by using He atom scattering (HAS) time of flight (TOF) spectroscopy 
to measure the dispersion of the low-energy adlayer-modes in combination with a theoretical 
analysis of the scattering data. This provides direct information on the force constants 
which couple the adsorbates to the substrate, the effective corrugation of the potential 
energy surface (PES) describing the motion of the particles on the substrate, and the force 
constants coupling the adjacent adsorbates. Although the first HAS-TOF measurements on 
rare gas overlayers were carried out almost twenty years ago, these early experiments were 
only able to detect the vertically polarized Einstein-like S-mode 0,0] • Only recently was an 
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additional mode found for Xe on Cu(llO) and shortly thereafter also for Xe on Cu(OOl) 



0, Cu(lll) and NaCl(OOl) ||TT| surfaces. This new mode was interpreted as being due 
to low-energy longitudinally polarized motion of the adsorbates. However, simple model 
calculations suggest a significant discrepancy between the fits based on this assignment 



TO] and empirical potential models used so far to model the adsorbate-adsorbate and 



adsorbate-substrate interactions. This has created some controversy in the literature ||12 
which calls for a critical reassessment of the HAS data on the low energy vibrations of 
physisorbed Xe monolayers. 

In this paper we present the analysis of the experimental data obtained for a commen- 
surate (v^ X ^/3)R30° Xe adlayer adsorbed on Cu(lll) and a similar system, viz. the 
hexagonal incommensurate Xe adlayer adsorbed on Cu(OOl). Measurements are reported 
on the dispersion curves and the excitation probabilities of two adlayer induced modes, the 
vertically polarized Einstein-like phonon branch over the entire surface Brillouin zone (SBZ) 
and the new in-the-surface plane polarized acoustic-like mode over the main part of the SBZ. 
In particular, for the latter mode in Xe/Cu(lll) it was possible to measure the zone-center 
phonon gap which provides particularly important information on the Xe-substrate potential 
energy surface, and to compare this result to the recently published data on the incommen- 
surate Xe-overlayer adsorbed on Cu(OOl) The first, tentative assignment of the modes 
observed in both systems was made by combining the symmetry selection rules applicable 
to the excitation of phonons of different polarization in the single phonon scattering regime 
with the results of lattice-dynamical analysis of the adlayers. This analysis showed that 
a good agreement between theoretical and experimental dispersion curves for Xe-induced 
modes could be obtained only with intralayer force constants derived from the Xe-Xe pair 
interaction unexpectedly softer than the one deduced from accurate gas phase potentials 
||r3| or bulk phonon data [|T^. As the use of this procedure was recently questioned |]T2[, fur- 
ther progress in understanding the HAS data for monolayers of Xe on Cu(lll) and Cu(OOl) 
necessitates resolving this controversy by carrying out a comprehensive and consistent anal- 
ysis of the measured spectral features. Here we demonstrate that additional support to 



the discussed assignments, and in particular of the dominantly longitudinal character of the 
observed acoustic mode, can be obtained by carefully examining the scattering intensities 
in the HAS TOF spectra. To this end we have performed detailed calculations of the ab- 
solute and relative excitation probabilities of the adlayer localized modes in the single and 
multiphonon scattering regime using a recently developed theoretical approach |T5|JT^ and 
compared them with experimental HAS TOF intensities. A good agreement of the results 
of calculations with the experimental data provides a strong and convincing argument in 
favor of the present assignments and therefore of the strong modification of the monolayer 
force constants. It thereby points to some peculiar characteristics of the Xe-Xe interaction 
in adlayers of monoatomic thickness which still await a microscopic interpretation through 
detailed calculations of the adlayer electronic and structural properties. 

The present article starts with a brief description of the experimental procedure and the 
preparation of Xe monolayers outlined in Sec. The experimental results are presented 
in Sec. |T| where a tentative assignment and interpretation of the modes observed in HAS 
from Xe monolayers on Cu(lll) and Cu(OOl) is discussed. The theoretical model used to 
interpret the HAS TOF intensities measured both in the single and multiphonon scattering 
regimes is described in Sec. The basic ingredients of the model relevant to the calcula- 
tions of HAS from adlayers are discussed and illustrated for the example of realistic He-Xe 
interaction potentials. In Sec. the experimental HAS TOF intensities are compared with 
the theoretical model predictions and physical implications of the good agreement are dis- 
cussed. Section ^ presents a way to determine the corrugation of the adlayer-substrate 
potential energy surface using the present approach. Finally, Sec. [VI I| summarizes all the 
relevant findings and the most important conclusions. Preliminary reports on some results 
of the present work have been published previously P,p!0[]. 
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II. EXPERIMENTAL 



The experiments have been carried out in a high resolution hehum atom scattering 
apparatus (base pressure 8 x 10^^^ mbar, fixed total scattering angle 6sd) described else- 
where [l^- The crystal was mounted on a home made crystal holder and the sample could 



be cooled down to 40 K with liquid helium. The sample temperatures were measured using 
a cromel-alumel thermocouple and could be stabilized to within ±0.5 K. The Cu(lll) crys- 
tal previously oriented to within 0.15° was cleaned by cycles of Ar ion sputtering (600 eV, 
1/iA/cm^) followed by annealing at 1100 K. Prior to the experiments described here the 
cleanliness and the structural quality of the surface were checked by XPS and LEED. The 
99.99% pure Xe gas was backfilled into the scattering chamber through a leak valve to a 
pressure of 4 x 10^^ mbar. During the preparation of the Xe monolayer structures the crys- 
tal temperature was kept at 70 K in order to avoid the formation of bi- and multilayers, 
which are stable only at temperatures below 68 K. For temperatures below 68 K the He- 
atom TOP spectra revealed the appearance of an additional loss, which could be assigned 
to vertically polarized vibrations of the second layer of Xe atoms. At the temperature of 
preparation of the Xe monolayer on Cu(lll) the commensurate-incommensurate transitions 



are not expected as they have been detected to occur at around 47 K |18|. A description of 



the preparation of the Xe monolayer on the Cu(OOl) surface has been given in Refs. 

III. EXPERIMENTAL RESULTS 

The right hand side of Fig. |l|.a shows the structure of the commensurate (v^ x -\/3)R30° 
monolayer of Xe atoms adsorbed on Cu(lll) surface [^,0] and indicates the two principal 



directions (azimuths) of the substrate crystal surface. The left hand side (LHS) shows the 
first SBZ of the substrate (dashed lines) and the two dimensional Brillouin zone of the 
adlayer (full lines). Figure ||.b shows an angular distribution of He atoms scattered from 
(V3 X V3)R30° Xe/Cu(lll) surface for incident wave vector h = 9.2 A-^ {Ei = 45 meV) 
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and the substrate temperature Ts=60 K along the [110]-azimuth relative to the substrate 
surface. The intensities are normalized to the specular peak height. In addition to the (1,0) 
diffraction peak, two additional, Xe (1/3,0) and (2/3,0) diffraction peaks of the order one- 



third were observed |T9[ . The sharpness of the peaks and relatively low background indicate 
the presence of a well-ordered, largely defect-free Xe-overlayer. 

Figure |^ shows an angular distribution along the substrate [100]-azimuth obtained by 
scattering He atoms from a monolayer of Xe atoms adsorbed on Cu(OOl) for incident wave 
vector ki = 5.25 A^^ {Ei = 14.36 meV) and Ts=52 K. This distribution indicates a well 
defined structure. The earlier LEED studies at > 77 K have identified a hexagonally 



ordered adlayer incommensurate with the underlying substrate p2| , p3| , so as that the [100] 
substrate direction lies along equivalent non-high symmetry directions for the two domains 
of the Xe monolayer, one rotated by 30*^ from the other. The geometrical structure of 



these domains is shown in Fig. 6 of Ref. |^ and Fig. 2 of Ref. [Q. The present out 



of-the-high-symmetry-plane measurements indicate sharp diffraction peaks commensurate 
with the orientation of the Xe overlayer along the Cu(OOl) [110] azimuth except for the two 
bumps at each side of the specular peak. Their presence is related to inelastic resonance 
processes involving the Cu(OOl) surface Rayleigh wave and intense nondispersive Xe mode 
with polarization perpendicular to the surface, as determined by the TOF measurements 
described below. These characteristics of He-angular distribution spectra from Xe/Cu(001) 
for Tg < 65K indicate a "floating" incommensurate Xe adlayer. Hence, in the case of both Cu 
substrates the Xe adlayers may be considered planar and periodic with hexagonal symmetry, 
irrespective of the (in)commensurability with the underlying substrate. The periodicity and 
symmetry of the adlayers is then reflected in their vibrational properties. The adlayer 
vibrational modes can be classified as dominantly in-plane polarized (longitudinal (L) and 
shear horizontal (SH)) and shear vertical (S) ||24|| . 

Figure ^ shows typical He atom TOF spectra for the (\/3 x v^)R30°Xe/Cu(lll) surface 
along the [112] substrate azimuth (i.e. along the TKxe direction of the superstructure), 
for 6sD = di + Of = 90.5° and three different He atom incident energies Ei spanning the 



transition from the single to the multi-quantum scattering regime. The TOF spectra have 
been converted from flight time to energy transfer scale. Arbitrary units are used for spectral 
intensities on the vertical axis. The spectrum at the lowest incident energy {Ei=9.9 meV) 
is typical of the single phonon scattering regime and is dominated by two well defined peaks 
at ±2.62 meV on the energy loss and gain sides of the TOF spectrum, respectively. Within 
the experimental error these energies do not change in the interval beyond AK = 0.1 
in the first SBZ of the superstructure. In accordance with previous works on noble 
gas atoms adsorbed on other substrates 0,^, this mode is assigned to the excitation of 
collective vibrations of Xe-atoms with a polarization vector vertical to the surface and is 
designated the S-mode. The lack of dispersion indicates that the frequency of the vertically 
polarized phonon is mainly determined by the adsorbate coupling to the substrate, with 
only a weak coupling to adjacent adsorbates. Deviation from a dispersionless behavior 
occurs only at the intersection with the substrate Rayleigh-mode [p5| , p6[ . The energy of this 
S-mode {huJs = 2.62 meV), is slightly larger than for the (llO)-face of Cu {huJs = 2.5 meV 
1^). The small, but significant deviation of 0.12 meV is consistent with a slightly deeper 
potential well for Xe on Cu(lll) and Cu(OOl) than on Cu(llO) 



In addition to the intense S-peaks the measured spectrum also reveals the presence of 
a weak but clearly resolved feature (labeled "L") near the elastic or zero energy loss line. 
The energy of this mode changes with the angle of incidence di and thus shows dispersion. 
The spectral intensity of this mode relative to that of the S-mode in each TOF measure- 
ment was found to decrease strongly with the magnitude of its wave vector so that for the 
He— s>Xe/Cu(lll) system the corresponding data points could only be obtained for parallel 
wavevector up to one third of the distance between the T and Kxe points in the first Bril- 
louin zone of the superstructure (interval shown in Fig. Since in the displayed TOF 
spectra the energy of the "L" mode is always significantly below that of the lowest surface 
phonon of the clean Cu(lll) surface |2^, this must be a pure Xe adlayer-induced mode 
which cannot couple to the substrate for wavevectors over a wide range of the SBZ. As 
for the (\/3 x ^/3)R30° Xe/Cu(lll) system the [112] direction has a high symmetry mirror 
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plane the vibrational modes are partitioned in two orthogonal classes |^0|. Two thirds of the 
modes are polarized in the sagittal plane (including the adlayer induced S- and L- modes). 
The remaining one third of the modes are polarized in the surface plane and normal to 
the mirror plane and designated shear horizontal or SH-modes. The three possible adlayer 
induced orthogonal modes with the wavevector in the [112] direction (c.f. Fig. |I].a) are thus 
characterized by either a combination of the components with S- and L-polarization or pure 
SH-polarization. Combining the symmetry selection rules pertinent to the probabilities of 
excitation of in-plane phonons at ideal surfaces (c.f. Refs. |3T1 , |32|| and Sec. |I^ below) with 



the fact that the data were recorded in the first SBZ of the superstructure and in the sagittal 
plane which coincides with the high symmetry plane of the Xe/Cu(lll) system, the obser- 
vation of the SH-mode under these experimental conditions can be ruled out. Hence, this 
mode is tentatively assigned to the longitudinal mode of the adlayer which is known to cou- 



ple to the scattered He atoms under similar conditions pTl32|. However, as demonstrated 



for NaCl, the SH-modes can be excited along a high symmetry direction in the second SBZ 



The other two spectra in Fig. ^demonstrate the transition from a single to a multiphonon 
scattering regime as Ei is increased. This transition takes place already at rather low He 
atom incident energies due to the very low excitation energies of the adlayer-induced S-modes 



whose vertical polarization gives rise to a strong projectile-phonon coupling (c.f. Sec. [rV| ). 
Although some single phonon features are still discernible at incident energy Ei = 21.4 meV, 
both spectra are dominated by a number of uniformly spaced peaks at energies ±n x 2.62 
meV. For Ei = 45.1 meV the true multiphonon scattering regime is reached because the 
intensity of the elastic peak is smaller than that of the multi-quantum S-peak for n = 2. 

Figure ^ shows three representative He-atom time-of-fiight spectra for the incommen- 
surate Xe monolayer on Cu(OOl) for three different He atom incident energies along the 
(100) substrate azimuth which lies halfway between the two high symmetry directions of 
the adlayer SBZ. In all essential aspects these spectra are similar to those shown in Fig. 2. 
They also exhibit strong dispersion of the "L" mode and the multiple excitation of S-modes 



at energies ±n x 2.71 meV. In addition, the Rayleigh mode (labeled RW) of the underlying 
Cu(OOl) substrate is also observed at low and intermediate energies Ei. The RW dispersion 
curves of clean Cu(lll) and Cu(OOl) surfaces are well known from the previous work p5| , p9| . 

It is noteworthy that for both adlayers the multiphonon lines are all, to within experi- 
mental error, located at integral multiples of a fundamental frequency us, 2.62 meV/h for 
Xe/Cu(lll) and 2.71 meV/h for Xe/Cu(001). At the first sight this seems to imply a very 
harmonic Xe-Cu potential since anharmonic shifts, which are expected to be negative, would 
produce overtone energies smaller than the corresponding multiples of the fundamental fre- 
quency uJs (c.f. Sec. [VI|) . However, the theoretical analyses of the Xe-metal interactions 
p7| , |28|j3^ show that the potential is notably anharmonic but, as we show in the next sections, 
the multiple spectral peaks can be explained by multiple excitation of delocalized phonon 
modes which involve the lowest harmonic states of many adatoms rather than a single higher 
anharmonic state localized on a single adatom. In this case there is no anharmonic shift as 
each multiphonon excitation is distributed over several Xe atoms in the overlayer. 

From up to about a couple of hundred of TOF spectra with different 6i the experimental 
dispersion curves were determined and are shown in Figs. ^ and ^. For both Xe/Cu(lll) 
and Xe/Cu(001) the vertically polarized S-mode exhibits negligible dispersion over the major 
part of the SBZ except at the point of avoided crossing with the substrate RW |3^. The most 



striking difference between the vibrational dynamics of the two adsorbate phases manifests 
itself in the dispersion of the "L" mode. The "L" mode for the commensurate Xe/Cu(lll) 
structure exhibits a zone center gap of about 0.4 ±0.1 meV whereas for the incommensurate 
phase the frequency at the zone center goes to zero linearly with the wave vector. 

To further aid the assignments of the modes in the He— >-Xe/Cu(lll) TOF spectra and 
analyze their dispersion we have carried out a full lattice dynamics calculation of the vi- 
brationally coupled (v^ x ^/3)R30° Xe/Cu(lll) system with Xe atoms placed in on-top 
sites on both sides of a 40 layer thick slab of substrate atoms and interadsorbate distance 
^Xe-Xe _ ^ j^2 A. The interaction between nearest-neighbor Cu atoms was accounted for by 
a single radial force constant (3'"^ = 28.0 N/m as obtained from a fit of the bulk Cu phonon 
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dispersion curves The other parameters describing the couphng of the Xe atoms to 

the nearest neighbor Cu substrate atoms was fitted to the dispersion curves, which yielded 
a radial force constant jS^'^ = 3.7 N/m and a tangential force constant a^'^ = 0.086 N/m. 
Assigning the longitudinal character to the observed "L" mode to comply with the above 
discussed symmetry selection rules, the interaction between the atoms in the adlayer could 
be described by a radial force constant f]^'^-^'^ = 0.5 N/m and a tangential force constant 



a 



Xe-Xe 



0. The results of the full calculation are shown in Fig. ^ in the next section 



and in Fig. ^ we present only the dispersion for the surface projected S- and longitudi- 
nal modes which reproduce the experimental data very satisfactorily. The radial Xe-Xe 
force constant fj^f^-^f^ = 0.5 N/m resulting from this procedure is, however, significantly 
smaller than the value predicted from the highly precise HFD-B2 gas-phase potential ||13|| , 



p^^fj^^ = 1.67 N/m, which produces a significantly steeper dispersion curve for longitudinal 
phonons denoted by dash-dotted curve in Fig. |^. 

In the case of incommensurate monolayer of Xe on Cu(OOl) it was only possible to set 
up a three-dimensional dynamical matrix describing the three vibrational modes localized 
in the adlayer, namely the S-, L- and SH- modes (c.f. Refs. |]7|,P^). Since the experimental 
data were taken halfway between the two high symmetry directions of the two dimensional 
hexagonal Xe adlayer Brillouin zone, the broken symmetry no longer forbids excitation of 
the SH-mode. However, since the calculated polarization vector of the SH-mode is nearly 
perpendicular and of the L-mode nearly parallel to the present azimuthal direction ||3^, the 



corresponding excitation probabilities of the SH-phonon are expected to be much smaller 
than those of the L-phonon (c.f. Refs. pl| , |32| and Sec. |rv| ). Hence, as in the case of the 
commensurate system, a longitudinal polarization is assigned to the observed low energy 
adlayer- induced acoustic "L" mode also in the incommensurate system Xe/Cu(001). The 
interadsorbate distance was fixed at d^^~^^ = 4.40 Aand the best-fit force constants in the 
[001] direction are: (3^^ = 3.8 N/m, /^^^-^^ = 0.42 N/m, and a^^'^^ = 0.012 N/m, which 
are similar to force constants for the Cu(lll) substrate. These results are also presented in 
Fig. P and reproduce the experimental data very well. For comparison the dispersion of the 
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L-mode calculated by using is also shown as a dash-dotted line and does not fit to 

the data. However, in the case of both Xe adlayers, the physical origin of the unexpected 
large softening of the radial Xe-Xe force constants introduced to reconcile the symmetry 
requirements with the experimental data remains unclear. A clue to this effect in the case of 
Cu(OOl) and Cu(lll) surfaces may be provided by their peculiar electronic structure which 
gives rise to surface electronic states with corresponding electronic wave functions extending 
far across the adsorbed Xe atoms . Alternatively, a delocalization of the electronic charge 



within the monolayer itself could give rise to softening of intralayer forces which then might 
explain the same phenomenon observed for an insulating substrate like NaCl(OOl) ||11|| . 

Recently, it was suggested that the observed "L"-mode could be interpreted as a SH- 
mode as this would be consistent with a thermodynamic analysis of the Xe/Cu(001) system 
This has created the controversy in the literature referred to in the Introduction and 



additionally motivated the present study. However, the model calculations of the HAS-TOF 
intensities reported in the next sections do not support such an interpretation. 

IV. THEORETICAL DESCRIPTION OF He ATOM SCATTERING FROM Xe 

ADLAYERS ON Cu(OOl) AND Cu(lll) 

Since the assignments of the peaks in the experimental TOF spectra discussed in the 
preceding section have been questioned ||12[ additional corroboration by theoretical argu- 
ments is called for. To this end we have carried out extensive calculations based on the 



recently developed fully quantum model of inelastic He atom-surface scattering |15|jl6|j39|j40 
which is especially suitable for scattering from adlayers Ill0| , |l9i^ . In this model, the single 
and multiphonon excitation processes can be treated on an equivalent footing without in- 
voking additional quasi-classical approximations for the scattered particle dynamics |^ , ^ 



A detailed description of the model was presented in Ref. [T^, so only its salient properties 
relevant to the calculations of HAS from adlayers are outlined here. 

In view of the HAS angular distributions characteristic of the present Xe/Cu(lll) and 
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Xe/Cu(001) surfaces (c.f. Figs. |l| and|^), the static corrugation of the He-surface interaction 
potential will be neglected so that the Xe overlayers are assumed to be flat and free of 
defects. The assumption of a planar and perfectly periodic Xe adlayer justifies the use of 
the dynamical matrix approach in the description of the vibrational properties of the surface 
in terms of phonon modes characterized by their parallel wave vector and branch index. The 
model Hamiltonian describing inelastic atom-surface collisions can then be cast in the form: 

H = Hr' + H^'' + V{r), (1) 

where 

= ^ + U{z) (2) 

is the Hamiltonian describing unperturbed motion of the projectile in the fiat static po- 
tential U{z) of the target. The projectile particle is characterized by its coordinate and 
momentum operators r = (p, 2;) and p = (P,^^), respectively, and mass M. Here p, z 
and P, Pz denote the parallel (lateral) and vertical (normal) to the surface components of 
r and p, respectively. H^^ is the Hamiltonian of the unperturbed phonon field which can 
be constructed once the dispersion and polarization of the vibrational modes of the Xe/Cu 
system are known, and V{y) is the dynamic projectile- surface interaction. The distorted 
waves (r I k) = (p, z \ K, A;^), which are the eigenstates of Hq""^^, are described by the quan- 
tum numbers K = P/^ and = Pz/^ denoting, respectively, the asymptotic parallel and 
normal projectile wave vectors far outside the range of the potential U{z). The correspond- 
ing unperturbed energy of the projectile is then given by Ek = h^iK^ + k1)/2M. Using 
the box normalization these eigenstates can be orthonormalized to satisfy (k | k') = 5k,k'5 
and are except for a phase factor equal to the unperturbed incoming and outgoing states 
satisfying the scattering boundary conditions [^. Explicitly, for a flat surface we have, in 
the coordinate representation, 
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where Ls and Lz are the quantization lengths in the parallel and perpendicular to the 
surface directions, respectively, and XkX^) satisfies the limit ^ oo) 2cos{kzZ + rj). 

A more detailed description of the distorted wave scattering formalism was given in Refs. 



3T|,|32|g|-|6| 



The theoretical description of the scattering event described by Hamiltonian is sought 
in terms of a scattering spectrum N{AE, AK), which gives the probability density that an 
amount of energy AE and parallel momentum hAK. are transferred from the He atom to 
the substrate phonon field. The particular choice of these two variables is dictated by the 
symmetry of the problem, the conservation of the total energy and, for a periodic surface, the 
conservation of the parallel momentum to within a reciprocal lattice vector G. Therefore, 
AE and AK completely determine the final state of the scattered particle provided its 
initial state ki = (Ki, kzi) is well specified. With these prerequisites the energy and parallel 



momentum resolved scattering spectrum is defined by |lT5| , |39| , |40| , p^ , |47| --^ 



N^,{^E, AK) = hm {^{t) I 5[AE - {H^ - £.)]5(^AK - P^'^) | v[/(t)), (4) 



where P^'' is the parallel momentum operator of the unperturbed phonon field, Si is the 
initial energy of the phonon field, AK = Ki — Kf, and | ^E'(t)) is the wave function of the 
entire interacting system. For calculational convenience, the phonon creation (annihilation) 
processes are assigned positive (negative) AE and AK since they refer to the phonon field 
quantum numbers. The corresponding experimental quantities have, however, opposite signs 
because they refer to the measured changes of energy and parallel momentum of the pro- 
jectile. The spectrum (^) is inherently normalized to unity and hence satisfies the optical 
theorem. It is also directly proportional to the experimental time-of-fiight (TOF) spectrum 
I6| . Based on the translational symmetry of the system we shall in the following describe 



the phonon modes propagating in Xe monolayers by their wave vector Q parallel to the 
surface, branch index j, polarization vector e(Q, j) and energy hujQj []7|,p^. The problem 
of localized modes characteristic of either mixed layers with broken or reduced translational 
symmetry |^ or of isolated adsorbates will be treated elsewhere. 
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Under the conditions in which the HAS experiments were carried out, the uncorrelated 



phonon exchange processes dominate over the correlated ones |]T5[. Then the angular resolved 



scattering spectrum is accurately represented by the expression: 

Ni.%^{^E, AK) = r ^I^e^K^^)--'^^^^)^! exp[2iy^^^(R, r) - 2iy^^^(0, 0)], (5) 

where r and R = (X, Y) are auxiliary variables occurring after the temporal and spatial 
Fourier representation of the energy and parallel momentum delta-functions, respectively, 
are introduced on the RHS of expression (^) [0, 2W^^^(Il, r) is the exponentiated Born ap- 
proximation (EBA) expression for the so-called scattering or driving function which contains 
all information on uncorrelated phonon exchange processes in the atom-surface scattering 
event. Its zero point value 2W^^'^{0, 0) = 2W^^^ gives the Debye- Waller exponent (DWE) 
in the EBA and the corresponding Debye- Waller factor (DWF), exp[— 2W^^^], gives the 
probability of the elastically scattered specular beam |5l|. Since in the EBA the correla- 
tions between two subsequent phonon scattering events are neglected, the expression for 
N^,^{AE , AK.) on the LHS of (|^) must be combined with the conservation laws for the 



total energy and parallel momentum |52 



The EBA expression for the scattering spectrum (^) holds irrespective of the form of 
the projectile-phonon coupling. However, it has been shown [^] that, for the projectile- 
phonon coupling to all orders in the lattice displacements, the higher order phonon exchange 
processes which involve only single phonon vertices and originate from linear coupling give 
much larger contribution to the scattering matrix than the non-linear many-phonon processes 
of the same multiplicity (c.f. Fig. 1 of Ref. and Fig. 1 of Ref. ||15|). Hence, in the 



present approach, only the linear projectile-phonon coupling will be retained, in which case 
the scattering function takes the form ||15|| : 



Here n{uQj) is the Bose- Einstein distribution of phonons of energy huQj at the substrate 
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temperature Ts and the symbols 

denote the on-the-energy and momentum-shell one-phonon emission (+) and absorption 
(— ) matrix elements or the probability amplitudes of inelastic He atom-surface scattering 
Tsj expressed through the corresponding off-the-energy-shell interaction matrix elements 



^fc^fc^f^^'^'^''' (s^^ below). In the present fully three dimensional calculations the wave 
vectors of real phonons exchanged in the collision will be restricted to the first SBZ of the 
superstructure (i.e. G=0) because the major part of the experimental data was recorded in 
this region. However, for integrated quantities such as the DWF, the relevant summations 
sometimes need to be extended beyond the first SBZ (see below). The matrix elements 
given by expression (|^) are normalized to unit particle current normal to the surface, jz = 
Vz/Lz = hkz/MLz. This can be easily verified if according to the box normalization the 
energy conserving (5-function in (|^ is converted into the Kronecker symbol following the 
identity 

2vr5(^KTQ,fci - EK,k, ± huQ^j) = ^X— SK,h.Q{kli), (8) 

Aiy Uz 

where k'^^ = ±2Ki ■ Q — + A^^j =F 2MujQj/h and 0(A;^j) is the step function restricting k^^ 
only to open scattering channels in which hkf /2M > 0. The factor of Lz appearing in (§) is 
canceled by the factor L^^ appearing in V'^^'^'''*^'^''-' which arises from the box normalization 
of the projectile wave functions (c.f. Eq. (|^)). This enables a straightforward summation 
over k'^ on the RHS of expression (^). The quantization lengths in the parallel directions 
also cancel out from expression (|^) after the final summation over Q. The Debye- Waller 
exponent is then obtained from Eq. (^) by a straightforward substitution R = and r = 0. 

In the following we shall adopt for V^'l^^^'^'''' the expressions obtained by taking the 
matrix elements of the vibrating part V"(r) of the total He-Xe/Cu potential Vtot(r) which 
will be modeled by a pairwise sum of atomic He-Xe interaction potentials f (r — ri) 



55| , |56[| . Hence: 
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V^(r) = (\4o*(r))„, = ^(r " ri - uj)) , (9) 

V 1 J vib 

where ri = (pi, zi) ranges over the equihbrium positions of Xe atoms and ui denotes the 
displacement of the 1-th Xe atom from equilibrium. The static part of Vjot(r), which is equal 
to the average of the pairwise sum in the brackets on the RHS of Eq. (|^), is then identified 
with the static atom-surface potential U{z) which is included in ifg"^* in Eq. (|^). The 
thus formulated scattering theory is equally appropriate to treat single and multiphonon 
scattering processes in HAS and incorporates the single phonon DWBA scattering theory 



reviewed in Refs. 31 and 32 clS cl special limit. 

To simplify the numerical calculations and in particular the treatment of the interac- 
tion matrix elements and the scattering function, we shall approximate U{z) by a suitably 
adjusted Morse potential: 

U{z) = /}(e-2"(--"o) _ 2e-"(^""°)), (10) 

where D, zq and a denote the potential well depth, position of the minimum and the inverse 
range, respectively. This approximation is justified in the range of energies of the present 
HAS experiments The values of these potential parameters appropriate to the two 



studied collision systems are given in Sec. ^ in which they are explicitly quoted as the input 
in the calculations of the scattering intensities for comparison with the experimental results. 



Following the findings of Ref . |58[ that the repulsive and attractive components of the pair 
potentials contribute to the vibrating part of the total dynamic atom-surface interaction, 
both components are included in the dynamic atom-adlayer interaction. Then the matrix 



elements for linear atom-phonon coupling acquire a simple form ||3T|j3^ : 



VjX"^'^'^' = u*(Q, j) ■ F(K' - K, K, h)SK',K^Q, (11) 

where u(Q,j) is the vector of quantized displacement of the adlayer atoms corresponding 
to the phonon mode (Q, j) defined below in Eq. (|15|). The parallel momentum conserving 
Kronecker symbols ^k'.k^^q arise as a result of the summation over adsorption sites of pair 
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potential contributions from all adsorbates in the periodic adlayer. The matrix element of 
the force F(K' — K, A;^, kz) exerted on the projectile by an atom in the adlayer is expressed 
as [0,0]: 



F(K' - K, A;;, h) = {XK I (^(K' - K), -d/dz)v{'K' - K, ^) | XfcJ, (12) 

where f (K' — K, z) is a two-dimensional Fourier transform of f (r) which appears in expres- 
sion (|l^) after taking the matrix element between the parallel components of the projectile 
wave functions (p | K') and (p | K). The scalar product appearing in expression (pUD 
gives the afore mentioned symmetry selection rules for excitation amplitudes of the various 
one-phonon scattering processes. In particular, for a purely SH-polarized phonon with a 
polarization vector in the surface plane and perpendicular to the mode wave vector Q in the 
first SBZ, the scalar product in (0) will be equal to zero, leading to a vanishing excitation 
probability. Expressions (|lT]) and (|12D also indicate that for small parallel momentum trans- 
fer to phonons the projectile-phonon coupling is strongest for vibrations with polarization 
vector perpendicular to the surface. 

The parallel or (K' — K)-dependence of f (K' — K, z) was modeled by introducing the 
cut-off wave vector Qc which gives an approximate upper bound on the parallel momen- 
tum transfer in the one-phonon exchange processes (Hoinkes-Armand effect [0]). Denoting 
K' - K = Q this yields: 

t;(Q, z) = Z}(e-M^-^o)e-«'/2Q? _ 2e-"(^-^o)e-QVQ?). (13) 

with the property ^(Q = 0,z) = U{z). The appearance of different effective cut-offs in 
the Q-dependence of the repulsive and attractive components of f (Q, z) in Eq. (pTSf ) (viz. 
V^Qc versus Qc) is a consequence of the different ranges of the repulsive and attractive 
components of the model pair potential v{r — ri) consistent with expression (|T0|) In the 
case of exponential potentials of range the value of Qc is approximately given by pl| : 



Qc = \-, (14) 



where Zt is the (energy dependent) average value of the He atom turning point in the surface 
potential U{z). Although the numerical evaluation of the matrix elements ( [T^ ) avoids the 
explicit introduction of Qc ||5^ the effect of the cut off in the parallel momentum transfer 



remains The matrix elements described by Eq. (|T^) were all calculated by using the 
pair potentials v{v — y\) to obtain f (Q, z), Eq. (|IB|), with the parameters corresponding to 
the two studied systems explicitly defined in the next section. Figures |^.a and |^.b illustrate 
the behavior of the one dimensional off-shell matrix elements (xfc?, | '^^(Q = 0, z) | Xfcz) 
and {xk'^ I —{d / dz)v{Cl = 0,2) | XfcJ which determine the parallel and perpendicular 
components of expression ([12D , respectively. The minima occurring in these plots are the 
consequence of the competition between the contributions coming from the repulsive and 
attractive components of f (Q = 0, 2;). 

As is clear from Eq. ([TT|), the calculations of the full interaction matrix elements also 
require the dynamical displacements ui of atoms in the adlayer. For Xe adsorption on 
Cu(OOl) the overlayer structure is incommensurate with the underlying substrate which 
necessitates an approximate but consistent description of the adlayer vibrational properties. 
Following the experimental evidence indicating a defect free and well structured monolayer 
(c.f. Fig. ^ we assume as in Sec. |ITT| a planar "floating adlayer" with perfect hexagonal 
symmetry in which the phonon dynamics takes into account only the averaged Xe-substrate 
potential perpendicular to the surface as well as inter-adsorbate interactions. In this case 
the polarization vectors and eigenfrequencies are solutions of a three-dimensional dynamical 
matrix ||2^. Hence, the adsorbate displacements can be expanded in terms of normal phonon 



modes of the adlayer in terms of phonon creation (Qq j) ^'^'^ annihilation (ciqj) operators: 

which also defines u(Q,j) = {fi/2MaNaUjQjy^'^e{Q, j) in Eq. (|llj). Here pj is the parallel 
component of the radius vector of the 1-th adsorbate in equilibrium, and Ma and Na are 
the mass and the number of Xe atoms in the adlayer, respectively. In the floating adlayer 
model the index j ranges over the S- (shear vertical), L- (longitudinal horizontal) and SH- 
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(shear horizontal) adlayer modes. By construction the floating adlayer phonon modes are 
locahzed in the adlayer and the polarizations of L- and SH-modes are strictly in the surface 
plane. The corresponding dispersion curves calculated with the adjusted force constants 
quoted in Sec. |T| are shown in Fig. ^. Consequently, the symmetry selection rules for 
probability of phonon excitation in a floating adlayer are more stringent than in the case of 
a dynamical matrix description which includes the substrate dynamics in some way. This is 
reflected in the calculated probability of excitation of an SH-phonon for in-the-sagittal plane 
inelastic scattering within the flrst SBZ of the superstructure. This probability turns out to 
be much smaller than for the L-mode for all azimuths, and is exactly zero along the two high 
symmetry directions TKxe and TMxe of the superstructure. These effects will be further 
discussed below in conjunction with the effects which the scattering potential imparts to the 
scattering intensities and in Sec. where we discuss a comparison between experimental 
and theoretical results. 

To calculate the required eigenfrequencies and the corresponding polarization vectors of 
the full dynamical matrix corresponding to the commensurate Xe/Cu(lll) system we have 
used the force constants given in the previous section. The results for the dispersion curves 
are shown in Fig. ^.a. Here we can also trace how each phonon mode of the composite system 
is localized at the surface (i.e. within the adlayer) and how typical surface modes may get 
delocalized for certain values of the wave vector. A measure of the adlayer localization of the 
S-, L- and SH-modes near the center of the flrst BZ of the superstructure is shown in Fig. 
^.b. Further important information concerning the ellipticity of polarization (perpendicular 
versus longitudinal) of the L-mode in the same region is illustrated in Fig. ^.c. 

The calculation of the force matrix elements, Eq. (0), can be performed analytically ||61 



using the wave functions of Eq. (|^) and the two dimensional Fourier transform of the pair 
interaction, Eq. (p^ ) ||58| . With the aid of these matrix elements and using Eqs. (0) and (||), 



we can obtain flrst order or distorted wave Born approximation (DWBA) transition proba- 
bilities for He atoms in Eq. (|^). Multiplying them by the corresponding Bose distributions 
we obtain the DWBA inelastic state-to-state reflection coefficients for one phonon scattering 
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3T| . Carrying out the summation over k'^ by making use of Eq. (H), fixing Osd = + 0f, 
Ei and Tg, and varying 6i we obtain tfie Q-dependent scattering intensities which can be 
directly related to the one phonon intensities in the TOF spectra (c.f. next section). Here 
we shall only illustrate some of their general and most interesting features. 

Figure ^ shows a plot of the DWBA scattering intensities for emission (energy loss) 
or absorption (energy gain) of a single dispersionless S-phonon in HAS from the "floating" 
Xe adlayer on Cu(OOl) surface. The calculated intensities reveal a relatively simple and 
expected structure as a function of the exchanged phonon momentum which is mainly due 
to simple properties of S-phonons characteristic of the "floating" adlayer model (absence 
of dispersion and the polarization vector localized in and perpendicular to the adlayer). 
On the other hand, due to the more complicated model of vibrational dynamics of the 
commensurate Xe monolayer adsorbed on Cu(lll) surface the analogous intensities shown 



in Fig. |TO|.b exhibit a more complicated structure despite the similarity in the corresponding 
potential parameters (for the magnitude of the latter see next section). Here the S-phonon 
polarization vector becomes delocalized over the first few layers of the Xe/Cu slab for the 
values of Q at which the S-phonon dispersion curve meets the dispersion curves of substrate 
phonons (c.f. Fig. ^.a). This makes the coupling of the He atom to S-phonons weaker in 
this region of the Q-space which is then reflected in the occurrence of pronounced minima 
in the scattering intensities. 

However, for nondispersive S-phonons the intensities of the measured single quantum 
loss or gain peaks acquire additional weight due to the multi-quantum interference between 
the emission of n ± 1 and annihilation of n nondispersive phonons of different wave vectors 
|6^ , |63[| . These additional contributions are automatically included in the EBA with account 



of recoil effects, as demonstrated in Ref. |6^. The difference between the single S-phonon 
loss and gain scattering intensities obtained in first order DWBA theory and in the EBA 



(which takes into account such interference processes) can also be visualized in Figs. [T^.a 
and |TD|.b for incommensurate and commensurate Xe adlayers, respectively. It is seen from 
the figures that for given projectile incoming energies these differences may be already quite 
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substantial. In particular, the noticeably larger EBA intensities for larger wave vectors are 
solely due to multi-quantum interference processes. This necessitates the use of the EBA 
in calculating the scattering intensities for comparison with the experimental data. On the 
other hand, we have found that the same type of renormalization of the scattering intensities 
of dispersive L-modes by emission and absorption of nondispersive S-modes gives a negligible 
effect. 

The differences in magnitude between the projectile transition probabilities in one phonon 
loss and gain processes can be best illustrated by repeating the above calculation for the 
scattering intensities but without multiplying the phonon loss and gain transition probabili- 
ties by the corresponding Bose distributions. The results of such a calculation for an S-mode 
emission and absorption are depicted in the insets in Figs, [l^.a and ^.b. These differences 



arise from and give a measure of the quantum recoil effects and are more pronounced at 



lower projectile incoming energies (c.f. Fig. (6) in Ref. [0). The trend that the results for 
loss processes give smaller values than for the gain processes reflects the fact that the phase 
space or density of states for transitions of the projectile to a state with lower energy (case of 
phonon emission) is smaller than to a state with higher energy (case of phonon absorption). 

Figure |rT|.a shows first order or DWBA intensities for HAS from L-phonons in the in- 
commensurate adlayer of Xe on Cu(OOl) surface for finite substrate temperature Tg which 
determines the temperature of the Bose-Einstein distributions for adlayer phonons in ex- 
pression (^. In the present "floating" adlayer model for the dynamical matrix the L-mode 
frequency follows acoustic dispersion ujq^i oc clQ for small Q (c.f. Fig. ||). The correspond- 
ing polarization vector e(Q, L) lies strictly in the adlayer plane and remains dominantly 
parallel to Q even outside the high symmetry directions of the BZ of the superstructure. 
Consequently, the L-phonon DWBA intensity for small Q (for which ujq^l Q) and finite 
Tg becomes proportional to the factor | Q ■ e(Q,L) p kTg/ujQi^ which saturates at a finite 
value for Q —>■ 0. As a result of that and the properties of the matrix elements depicted 
in Fig. |^.a the L-phonon scattering intensity in the case of the incommensurate monolayer 
exhibits a maximum near the zone center. The minima in the intensity occur for those val- 
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ues of Q at which the on-the-energy-shell counterparts of these matrix elements go through 
a minimum or zero. All these features clearly manifest themselves in the Q-dependence of 
the L-phonon scattering intensity shown in Fig. |TT].a. In contrast to the L-mode behavior, 
the polarization of the SH-mode remains dominantly perpendicular to Q also outside the 
high symmetry directions of the first SBZ of the superstructure Therefore for ideally 



structured "floating" adlayers it will always give a much smaller contribution to the intensity 
of the one phonon processes for HAS in the sagittal plane and Q restricted to the first SBZ 
(c.f. next section). 

Figure |TT|.b shows the single L-phonon HAS intensity for Xe/Cu(lll) system as a function 
of the phonon wave vector on the same scale as in Fig. [l^.a. Apart from the trends leading 
to zero scattering intensities for some isolated Q-values due to the behavior of the off- 
shell matrix elements shown in Fig. ^.a, which are common to both incommensurate and 
commensurate Xe layers, some basic differences with respect to the Xe/Cu(001) system 
can be observed. The commensurability of the adlayer with the substrate gives rise to 
non-vanishing Xe-Cu shear stress force constants entering the full dynamical matrix of the 
Xe/Cu slab, thereby producing two important effects regarding the phonon dynamics. First, 
it causes the appearance of a zone center phonon gap in the dispersion curves of L-phonons 
and hence the intensity factor | Q-e(Q, L) p kTg/ujQ j^ oc for small Q because co'q=o,l 7^ 0. 
Second, the polarization vector of the L-phonons is no longer constrained to the surface plane 
but for the values of Q at which the L-mode and substrate modes are degenerate it also 
acquires a component in the direction perpendicular to the surface and its localization to the 
adlayer is reduced (c.f. Figs. |^.b and ^.c). For a completely in-surface-plane polarization 
of the L-mode and finite Tg the first effect causes a drop in the scattering intensity which 
is quadratic in Q towards the zone center. The beginning of this trend is clearly visible 
in Fig. |TT].b. However, by increasing Q the second effect begins to play a role and, since 
the coupling to perpendicular vibrations is much stronger than to the parallel ones, the L- 
phonon scattering intensity rises rapidly up and reaches a maximum in two spikes near the 
zone center. Hence, the interplay between the parallel and perpendicular polarizations, or the 
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ellipticity of the L-mode polarization in the commensurate Xe/Cu(lll) system, introduces 
fast variations of the scattering intensity as a function of the exchanged phonon momentum. 
Regarding the role of the SH-modes of the slab on the HAS intensities, we find that the 
situation here is completely analogous to the case of a fioating Xe adlayer, i.e. within the 
present model description their coupling to He atom is generally negligible for in-sagittal- 
plane collision geometry and Q lying in the first SBZ, and is strictly zero along the high 
symmetry directions. However, this conclusion does not apply to the magnitude of the 
Debye- Waller factor because the Debye- Waller exponent is obtained by carrying out the 
momentum transfer summations over the entire first SBZ and also beyond if the coupling 



matrix elements are strong there (c.f. Figs. [T^ and |TT]). In fact, for the matrix elements of 
appreciable magnitude in the second SBZ the SH-mode can produce even larger contributions 
to the DW exponent than the L-mode due to its lower excitation frequency. This situation 



is illustrated in Fig. |2. 

An important point to be observed in connection with Figs. ^and|Tl|is a relatively large 
difference between the one phonon loss and gain scattering probabilities at lower projectile 
incident energies. This is due to an interplay between the projectile recoil and temperature 
effects. The amplitude for a transition of the projectile to a state with higher energy is 
larger than to a state with lower energy (c.f. insets in Figs. |l^.a and [l^.b) because of the 
larger phase space for the states of higher energy. Thus, for the same incoming energy the 
matrix elements | V(— ) p for phonon absorption in expression for the scattering function 
will be generally larger than | V(+) ^ which describe phonon emission (c.f. also Fig. 6 in 
Ref. [|19[)- On the other hand, the temperature effects entering expression (H) through the 
corresponding Bose- Einstein distribution factors n and {n + 1) for phonon absorption and 
emission, respectively, act just in the opposite direction because n < (n -|- 1) . Hence, the 
total scattering intensity for phonon emission or absorption calculated from the scattering 
function (10) depends on the trade off between these two effects. 
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V. COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 



A real test of the assumptions underlying the present model interpretation of HAS from 
Xe monolayers on Cu(OOl) and Cu(lll) surfaces should come through a comparison of rela- 
tive excitation intensities of the various modes in the experimental and calculated scattering 
spectra. This is equally relevant for the single and multiphonon scattering regimes. 

In this section we first calculate the relative intensities of the various adlayer modes in the 
single phonon scattering regime using the approaches described in the preceding sections. 
The parameters characterizing the He-Xe potential in expression (|r^) are given by: 

D = 6.60 meV, a'^ = 0.8202 A, and zq = 3.49 A, 
for the commensurate system Xe/Cu(lll), and 

D = 6.40 meV, a'^ = 1.032 A, and Zq = 3.6 A, 
for the incommensurate system Xe/Cu(001). These are not the parameters obtained from 
pairwise summation of pure He-Xe gas-phase potentials which yields: 

Dgas = 7.2 meV, a~^^ = 0.77 A, and Zo^gas = 3.51 A, 
but slightly modified ones to produce a softer He-surface potential. The present set of 
parameters for the Xe/Cu(lll) system also differs from the one quoted in Ref. |1^ which 
is a result of additional consistency requirement imposed to obtain a unified set describing 
the TOF spectral intensities equally well in the single and multiphonon scattering regimes. 
Such a necessity for modification of the sum of gas-phase pair potentials is not uncommon 



in HAS studies (c.f. Refs. ||60|j65f| ) and here it is also necessary to consistently reproduce the 
relative peak intensities in the measured TOF spectra. 

A comparison of the experimental and calculated spectral intensities of S- and L-modes 
in the single phonon scattering regime of HAS from Xe/Cu(lll) is shown in Fig. |13|. Here 
we note that our calculations always yield smaller intensities for the elastic peaks compared 
with the experiment as the latter also includes contributions from diffuse elastic scattering 
from defects not accounted for by the present model. Hence, this missing component has 
been added to the elastic peak intensity because in combination with the finite peak width 
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it can also contribute to the background intensity of the L-peaks. With this proviso a good 
agreement between experimental and theoretical results is achieved which illustrates the 
consistency of the present interpretation of the inelastic peaks in the TOF spectra. 

Since in the floating adlayer model applied to the Xe/Cu(001) system the interference 
between the adlayer modes and the substrate Rayleigh wave (RW) cannot be obtained, 
we have selected the TOF spectra in which the peaks assigned to S- and L-modes could 
be maximally separated from those of the RW and then calculated only the adlayer mode 
intensities in the EBA by neglecting the S-mode frequency shifts and delocalization occurring 
at avoided crossing with the dispersion curve of the RW. Along the direction of measurement 
the polarization of the SH mode is not strictly perpendicular to its wave vector, and this in 
principle could give rise also to SH-mode-induced structures in the TOF spectra. However, 
in this direction the present theoretical analysis gives the SH-mode intensity of the order 
of only five percent relative to the contribution from the L-mode. Only the introduction of 
coupling of the SH mode to the modes of the underlying substrate could induce a breakdown 
of symmetry leading to a removal of such stringent selection rules and thereby to a larger 
SH-mode excitation probability. The calculated relative intensities of the S- and L-modes 
for the Xe/Cu(001) system in the single phonon scattering regime, but with inclusion of 
multi-quantum interference of nondispersive S-modes, are compared with the experimental 
data in Fig. |l^ Given all the approximations used in the calculation, it is seen that the 
model reproduces the relative TOF intensities quite satisfactory. The only exception occurs 
on the energy loss side in the lower panel of Fig. |l^ where the afore mentioned avoided 
crossing between the substrate RW (energy loss at 3.05 meV) and the adlayer S-mode takes 



place, with the effect of S-mode frequency shift and intensity enhancement. In Fig. |4 
we also demonstrate how the use of the unsoftened force constants derived from the three- 
dimensional Xe-Xe gas phase potential |T^ produces the position and intensity of L-peaks 
for which the disagreement with the experimental data is evident. On the other hand, 
the SH-peaks calculated for the present scattering geometry by using the same gas phase 
potentials are of negligible relative intensity to be experimentally observable although their 

26 



frequency may coincide with that of the measured acoustic mode. Hence, with the present 
mode assignments and the corresponding model dispersion relations based on the softened 
radial intralayer force constants we can consistently describe the HAS-TOF intensities for 
the Xe/Cu(001) system as well. 

As the coupling of He atoms to S-modes is much stronger than to the L-modes, as il- 



lustrated in Figs. [T^ and |Il|, the multiphonon scattering spectra will be dominated by a 
series of multi-quantum S-peaks. All other dispersive modes may only add weak structures 
on top of this basic one. Eventually, these structures will turn into a broad Gaussian-like 
background (cf. Eq. (73) in Ref. [|15| and Eq. (9) and Fig. 3 in Ref. [0) in the limit 
of high incident projectile energies. The multiphonon scattering spectra from Xe/Cu(lll) 



and Xe/Cu(001) adlayers have been studied in detail in Ref. [0] and for the sake of com- 
pleteness we here briefly illustrate their features in Figs. |l^.a and |T5|.b, respectively, by 
consistently employing the potentials from the corresponding single phonon calculations. 
Figure ^.b is interesting in that it illustrates the behavior of the experimental scattering 
spectra of He— s>Xe/Cu(001) collisions in the single phonon scattering regime regarding the 
modes which weakly couple to He atoms. This is shown by the appearance of a Rayleigh 
wave-induced hump and an L-mode-induced shoulder near the elastic line. However, this 
scattering regime simultaneously appears to be a multiphonon one for the S-modes whose 
coupling to the scattering He atoms is much stronger. Concerning the agreement between 
experimental and theoretical results we observe that although the multiphonon TOF spectra 
of the incommensurate Xe/Cu(001) system can be relatively well reproduced, except for the 
elastic line which was not corrected for diffuse scattering contribution in the theoretical plot, 
the agreement for the commensurate Xe/Cu(lll) system is better. Interestingly enough, in 



the latter case (Fig. |T5[a) the diffuse elastic scattering correction for the elastic line was 
unnecessary due to the true multiphonon character of the spectrum at this higher incident 
energy. This spectrum can be viewed as a convolution of a series of well defined equidis- 
tant peaks, signifying the uncorrelated multiple emission and absorption of nondispersive 
S-phonons and not the overtones (for overtone frequencies see next section), with a broad 
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background arising from the multiple excitation of L- and SH-phonons (whose polarizations 
in this regime are no more constrained by the symmetry selection rules), and also from the 
substrate surface projected modes which may couple to the scattered He atoms. As implied 
in expression (^), the phonon absorption processes of any kind can take place only for Tg > 0. 

It is also noteworthy that the resulting multi-quantum S-mode intensities in the spectra 
shown in Figs. |l5[a and |l^.b do not generally follow the Poisson distribution. This arises as 
a consequence of the projectile recoil effects and the dependence of the magnitude of each 
peak maximum on the exchanged parallel momentum, which both act so as to prevent the 
appearance of a simple Poissonian structure. These features were discussed in detail in Refs. 



VI. CORRUGATION OF Xe-Cu(lll) POTENTIAL ENERGY SURFACE 

The value of the zone center gap for the longitudinal phonon mode in the case of Xe 
adsorbed on Cu(lll) surface also provides information on the Xe/Cu(lll) potential en- 
ergy surface. Since, in the past, precise information on the lateral corrugation for noble 
gases adsorbed on metal surfaces has been missing, an approach based on the work by 
Steele has frequently been used in applications [^] in which the periodic adsorbate- 



substrate potential K(r) is obtained by a summation of Lennard- Jones pair potentials 
V{r) = 4e[(cr/r)^^ — (a/r)^] describing interaction between noble gas atom and substrate 
atoms. The resulting sum can be expressed in the form 

V,ip, z) = e[Vo{z) + / V^[z) X: cos(G, ■ p)]. (16) 

1=1 

where now p and z denote the adsorbate parallel and perpendicular coordinate, respectively, 
the latter being measured from the first layer of substrate atoms. The symbols ±Gj, i = 1 
to 3, denote the six shortest reciprocal lattice vectors, and the definition of Vq and Vi is given 
in Ref. |^^. The original expression (|16|) was derived in Ref. for the case / = 1 and 



the extra factor / which may differ from unity has been added following Ref. [B3] to adjust 
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the corrugation without changing the adsorption energy. Since this type of the potential 
has been used in molecular dynamics simulations of sliding monolayers we can estimate 
its relevance to the present problem even without invoking the precise information on the 
adsorption sites. Thus, setting a to 3.487 A, the arithmetic mean of axe-Xe (describing 
the Xe-Xe interaction), the Cu nearest neighbor distance pU| and putting e = 19 meV, 



yields a Xe-surface potential with a total depth of around 200 meV, a reasonable value P7 
For the interaction between noble gas atoms and metal surfaces with highly delocalized 
conduction electrons the parameter / has previously been reduced to values of the order 0.1 
to account for small corrugations observed experimentally, and values as small as 0.03 have 
been proposed for the case of Kr on Au(lll) 



Although we have not used the parametrized expression ([T6| ) in our dynamical matrix 
analysis of Xe adlayer vibrations, we would like to point out that for the potential parameters 
appropriate to Xe on Cu(lll) (see above) a value of / = 0.1 yields a zone center phonon gap 
energy of 0.25 meV, significantly smaller than the value reported here. In addition, using 
this potential energy surface (PES) the computed S-phonon energy amounts to 4.55 meV 
(i.e. much larger than the present experimental value hujs = 2.62 meV), and the anharmonic 
shift of the fifth overtone to ~0.5 meV. The same conclusion can be also derived by using 
the Xe-Cu(lll) potentials parametrized in Refs. |p7| , |28[| . To get the correct curvature of 



the PES at the adsorption site one should either change by a substantial amount the depth 
of the potential ([T6|) in the normal direction (governed by e), or modify its width which is 
controlled by cr, or both. Drastic changes in the potential depth leading to the experimental 
vibrational energies are physically implausible in view of the correct fitting of the desorption 
energies. On the other hand, the correct curvature can be achieved by a modest, 10% 
softening of the repulsive component of the potential. This is consistent with the softening 
of the L-mode force constants discussed in Sec. ^ and the modification (also in the direction 
of softening) of the He-Xe adlayer potential noted in Sec. |V|. All three features, which derive 
from independent analyses of the experimental data, unambiguously point towards the effect 
of softening of pair interactions involving Xe atoms in adsorbed monolayer phase. 
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In the previous report of the present results for Xe/Cu(lll) [|T0[ an incorrect reference 
was made to the work of Cieplak, Smith and Robbins in which the resuhs of molecular 
dynamics simulations for Kr adsorbed on Au(lll), i.e. not Xe on Ag(lll), were reported. 
Considering the analysis presented in the last paragraph, we believe that the PES used in 
that work [^] may be too weakly corrugated when the factor / is reduced to below the value 
of 0.1. Since the molecular dynamics results on friction of sliding noble gas adlayers show 



a pronounced dependence on the strength of lateral corrugation l6gJ7y] it would be highly 
desirable to experimentally determine the zone center gap of the longitudinally polarized 
phonon modes of Kr and Xe adsorbed on Ag(lll) or Au(lll). Unfortunately, however, 
such measurements are hampered by the fact that on these substrates both noble gases 
form incommensurate overlayers in which the zone center phonon gap is zero because of the 
translational invariance. 

VII. SUMMARY AND CONCLUSIONS 

In this work we have carried out a comprehensive comparative experimental and theo- 
retical study of low energy dynamics of monolayers Xe on Cu(lll) and Cu(OOl) surfaces by 
utilizing the HAS-TOF spectroscopy and a recently developed quantum theory of inelastic 
HAS from surfaces which treats single and multiphonon scattering processes on an equiv- 
alent footing. The inelastic HAS data were obtained for a wide range of initial scattering 
conditions {Ei,6i,Ts) spanning the single and multiphonon scattering regimes. The data 
have been carefully analyzed by combining the dynamical matrix approach for description 
of the adlayer vibrational dynamics with the developed scattering theory. 

For both substrates the angular distributions (diffraction spectra) of scattered He atoms 
signified hexagonally ordered Xe monolayers, the commensurate (v^ x ^/3)R30° adlayer in 
the case of the Cu(lll) surface p0|,|21|] and the incommensurate one in the case of the Cu(OOl) 



surface [22,23[. The measured phonon dispersion relations deduced from the TOF spectra 



for the two types of monolayers exhibited great similarities in the case of the nondispersive 
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optical-like S-phonons localized in the adlayer. However, a striking difference occurred in the 
case of the very soft dispersive phonon branch labeled " L" . In the incommensurate phase it 
is genuinely acoustic-like whereas in the commensurate phase it exhibits a frequency band 
gap at the zone center. 

The tentative mode assignment was first made by employing the symmetry selection 
rules applicable to one-phonon excitation processes in HAS in the sagittal plane in combi- 
nation with the dynamical matrix description of the adlayer vibrational modes. Since in the 
case of the Xe/Cu(001) system no signature of a possible high order commensurate phase 
has been found, the dynamical matrix treatment of the two adlayers was markedly different. 
On Cu(OOl) it was only possible to treat incommensurate Xe adlayer as an ideally struc- 
tured hexagonal monolayer floating on a structureless substrate surface with adlayer modes 
completely decoupled from those of the substrate. On the other hand, for the monolayer 
of Xe on Cu(lll) it was possible to construct a full dynamical matrix of the vibrationally 
coupled Xe/Cu(lll) system. Despite these differences it was possible to establish a unified 
model interpretation of the observed modes in both systems in terms of a consistent set 
of adjusted Xe-Xe and Xe-Cu force constants. In this picture the S-mode is dominantly 
adlayer localized and FT^ or vertically polarized. In order to reconcile the symmetry se- 
lection rules for excitation of in-plane-polarized adlayer-induced modes in HAS from ideally 
structured adlayers with the dynamical matrix description of the modes, we had to assign 
a dominantly longitudinal polarization to the "L"-mode (L-phonon), with a possibility of 
a relatively strong vertical or 2;-admixture (up to ~ 15%) near the avoided crossings with 
other surface projected modes. By consistently carrying out this procedure we found that 
the radial Xe-Xe force constants determining the dispersion of the thus assigned longitu- 
dinal modes turned out much softer than by directly applying the sophisticated HFD-B2 
Xe-Xe gas phase pair potential ||13|. Such a modification of the interadsorbate interactions 
still awaits its interpretation through a detailed calculation of the electronic structure of 
Xe adlayers on Cu(lll) and Cu(OOl). The magnitude of the band gap in the dispersion of 
L-phonons in the commensurate Xe/Cu(lll) phase was used to obtain information on the 
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corrugation of the Xe atom-substrate potential in this system . 

The assignment of the observed acoustic modes to shear horizontal or SH-modes in the 
first SBZ, which were the subject of recent controversy [|1^, was ruled out for both systems 
first by invoking the symmetry selection rules (c.f. Sec. 0) for the excitation of the various 
in-plane modes in HAS. Following these arguments, in the present model in which the ideal 
structure of Xe adlayers is an essential ingredient, the SH-mode can give only a negligible 
contribution to the scattering intensity relative to the longitudinal mode. 

These interpretations were then corroborated by detailed theoretical analyses of the scat- 
tering intensities in the HAS-TOF spectra. We have first analyzed and demonstrated the 
difference between the interaction matrix elements coupling the scattering He atoms to the 
perpendicular and parallel to the surface vibrations of the adlayer for small parallel momen- 
tum transfer and variable incoming projectile energy. Combining this with the calculated 
properties of the polarization vectors of the adlayer modes, we were able to obtain the one- 
phonon HAS intensities of the adlayer modes and compare them with experiments. A good 
and consistent agreement between experimental and model theoretical results for the S- and 
L-modes was obtained in all aspects of the measured data for both types of Xe adlayers. 

The theory was then extended to calculations of the scattering spectra in the transition 
from the single to the truly multiphonon scattering regime which is characterized by the value 
of the corresponding Debye- Waller exponent larger than unity. Interestingly enough, for the 
S-modes, which strongly couple to the scattered He atoms, the multiphonon regime is already 
reached for the scattering conditions which still favor the single L-phonon scattering. The 
interplay of these two types of couplings gives rise to the characteristic spectral shapes of the 
multiphonon HAS TOF spectra. A very good agreement between the calculated results and 
the TOF data in this regime gives further support to our earlier assignments, the consistency 
of the model description of the low energy dynamics of the two distinct monolayer phases of 
Xe, and the adequacy of the present treatment of single and multiphonon He atom scattering 
from these adlayers. 
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FIGURES 

FIG. 1. (a) Right panel: Structure of (-\/3 x -\/3)R30° monolayer of Xe atoms (shaded circles) 
on (111) surface of Cu crystal, with two high symmetry directions (azimuths) in the substrate 
surface plane denoted. Left panel: Two dimensional Brillouin zones of the Cu(lll) surface (dashed 
lines) and of Xe adlayer (full lines), (b) He-atom angular distribution along the [110] azimuth of 
Cu(lll) from Xe (\/3 x \/3)R30° monolayer for incident wave vector ki = 9.2 A-^ (-Ej=45 meV) 
and surface temperature 60 K. Peaks normalized to the height of the specular peak. 

FIG. 2. He-atom angular distribution along the [100] azimuth of the substrate from an incom- 
mensurate hexagonal monolayer of Xe atoms adsorbed on Cu(OOl) for ki = 5.25 {Ei = 14.36 

meV) and substrate temperature 52 K. 

FIG. 3. Series of measured HAS TOF spectra for a monolayer of Xc on Cu(lll) along the 
[112] direction of the substrate surface for three representative He atom incident energies. 6i and 
9sD denote the angle of incidence and the fixed total scattering angle, respectively. Vertical scale 
measures relative peak heights in arbitrary units. Scattering parameters shown in the insets. 

FIG. 4. Scries of measured HAS TOF spectra for a monolayer of Xc on Cu(OOl) along the 
[100] direction of Cu surface for three representative He atom incident energies. Vertical scale in 
arbitrary units. Scattering parameters shown in the insets. 

FIG. 5. Phonon dispersions for Xe/Cu(lll) surface along the [112] direction relative to the 
substrate as determined by HAS (full circles). The solid line denotes the best fit achieved for the 
longitudinal (L) mode in the Xe adlayer and the dashed-dotted line is the result for the L-mode 
using the gas phase Xe-Xe potential. The theoretical dispersion curve for the vertically polarized 
S-mode is marked by the long dashed line and of the Rayleigh phonon and the projected bulk 
phonon edge of the Cu substrate by the full thin and dotted line, respectively. For force constants 
see main text. 
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FIG. 6. Phonon dispersions for Xe/Cu(001) surface along the [100] direction relative to the Cu 
substrate as determined by HAS (full circles) . The full curve represents the best fit achieved for the 
longitudinal (L) mode, the flat dashed curve shows a theoretical dispersion curve for the S-mode, 
the dot-dashed line indicates the position of the L-mode using the gas phase force constant. Dashed 
curve indicates the position of the Rayleigh wave on the Cu(OOl) surface. For force constants see 
main text. 

FIG. 7. (a) Off-shell matrix elements of the potential which couples the projectile to in-plane 
surface vibrations (see text) for Q = and three different values of the incoming perpendicular 
wave vector: k^i = 3.5 A(full line), kzi = 7.0 A(dashed line) and kzi = 14.0 A(dotted line), as 
function of scattered wave vector k'^. Potential parameters corresponding to Xe/Cu(lll) system, 
(b) Same for projectile coupling to perpendicular vibrations. 

FIG. 8. Calculated dispersion curves of S-, L- and SH-phonons (from top to bottom) over 1/6 
of the two dimensional first Brillouin zone, i.e. between two equivalent high symmetry directions, 
corresponding to a floating Xe adlayer on Cu(OOl). For force constants parametrization see main 
text. The angle (j) is measured relative to TMxe direction of the Brillouin zone of the superstructure. 

FIG. 9. (a) Dispersion curves for Xe/Cu(lll) system along the [112] direction calculated using 
the dynamical matrix approach (see text) with Xe atoms on-top of Cu atoms. Force constants 
parametrization same as in Ref. ||l^. Note the S-, L- and SH-mode dispersion curves detached 
from the bulk quasi-continuum. (b) Surface localization of Xe induced S-, L- and SH-modes on 
the adlayer expressed through the sum of the components of the respective polarization vectors in 
the adlayer. Numbers below symbols A denote frequencies (in meV) of the S-mode. (c) Ellipticity 
(vertical vs. longitudinal polarization) of the L-mode in Xe adlayer on Cu(lll) substrate as function 
of phonon wave vector. Percentage above full squares denotes surface localization of the L-mode. 
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FIG. 10. (a) Calculated DWBA and EBA scattering intensities for emitting one S-phonon 
(lower and upper full curves, respectively, marked by (+)) and absorbing one S-phonon (lower 
and upper dashed curves, respectively, marked by (— )) in HAS from Xe/Cu(001) as function of 
exchanged phonon wave vector Q for the scattering conditions as denoted. For potential parameters 
see main text. Inset: Corresponding projectile transition probabilities in one phonon loss and gain 
processes denoted by full and dashed lines, respectively, (b) Same for Xe/Cu(lll) system. 

FIG. 11. (a) Calculated DWBA HAS intensity for L-phonon emission (full curve) and L-phonon 
absorption (dashed curve) in Xe/Cu(001) as function of exchanged phonon momentum. For inter- 
pretation of the maxima and minima see main text, (b) Same for Xe/Cu(lll) system. 

FIG. 12. Plots of the various contributions to the Debye-Waller factor pertinent to 
He— >Xe/Cu(001) collision system as function of substrate temperature Tg and for the scatter- 
ing conditions as denoted. Dotted line: only S-phonons included; dashed line: S- and L-phonons 
included; full line: S-, L- and SH-phonons included. Note the logarithmic scale on the ordinate 
axis. Here the Q summation in Eq. (^) to obtain the DW exponent was carried out over the first 
and second surface Brillouin zone of the superstructure. 

FIG. 13. Comparison of experimental and calculated S- and L-mode intensities in the sin- 
gle-phonon scattering spectra typical of He— >Xe/Cu(lll) collisions. 

FIG. 14. Comparison of experimental and calculated S- and L-mode intensities in the sin- 
gle-phonon scattering spectra typical of He— >Xe/Cu(001) collisions. Expanded contours in the 
topmost panel show the positions and relative intensities of the L-mode peaks calculated with 
the present softened Xe-Xe radial force constants (full curve) and unmodified gas phase poten- 



tial-derived 1 13 1 force constants (diamonds). Similar differences are obtained for other two spectra. 



SH-derived peaks not discernible on the present scale. 
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FIG. 15. (a) Comparison of experimental and calculated EBA multiphonon scattering spectra 
typical of He^Xe/Cu(lll) collisions using the same potentials and dispersion relations as in Fig. 
|l3|. (b) Same for He— >Xe/Cu(001) collisions using the potentials and dispersion relations as in Fig. 
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